Point Charge Dynamics Near a Grounded Conducting Plane 
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The classic image problem in electromagnetism involves a grounded infinite conducting plane and 
a point charge. The force of attraction between the point charge and the plane is identified using an 
equivalent-field picture of an image charge with opposite sign equidistant behind the plane resulting 
in a force of attraction between the original charge and the plane. If the point charge is released 
from rest it will reach the plane in a time r. This time r has not been calculated correctly up to 
now. Clarification of the inconsistency is presented along with a correct solution to the classic image 
problem. Other electromagnetism problems are mentioned with attractive l/r" forces (where n > 
1). Such situations arise between a point charge and a line charge (l/r), between a line charge and 
a point dipole (1/r^), between a point charge and a point dipole (l/r''), and between a point dipole 
and a second point dipole (l/r'*). 



I. INTRODUCTION 

Point masses under the influence of Newtonian gravity 
and point charges under the influence of Coulombic at- 
traction experience a l/r ^ force. Particles in this picture 
will approach one another according to Newtonian (non- 
relativistic) dynamics and if released from rest will reach 
the force center in a time that is typically shown to be 
proportional to the initial separation to the three-halves 
power. This type of solution appears in various places in 
textbooks as an application of gravitational forces [l|, 
and in electromagnetism textbooks surrounding the clas- 
sic image problem0-Q. Unfortunately, the solution lead- 
ing to this result is plagued with superluminal (i.e. faster 
than light speed) transport and hence cannot be formally 
correct. 

In addition to rclativistic effects, action-at-a-distance 
type forces require propagation of information about the 
force. News of changing particle positions must prop- 
agate to the other particle before any changes can be 
realized and the news travels at a fixed (finite) speed, 
namely, the speed of light. Care must be taken to prop- 
erly account for these retardation effects. The ultimate 
effect is to extend the time it takes for the particle to 
reach the force center as compared with an assumption 
of instantaneous updating. 

The electromagnetic category of problem is discussed 
in this paper, but the same formalism can be applied to 
the gravitational problem as well. The solution discussed 
below for problems of this type incorporates Einstein's 
special relativity and includes the effects of field retarda- 
tion. In Section [n] the classic image problem is discussed 
and solved taking special relativity into account. After 
all, when particle separation approaches zero with an in- 
verse squared force law, infinite kinetic energies are possi- 
ble. This has all the earmarks of superluminal dynamics. 
Consistent results without any superluminal dynamics 



are presented in Section (|TI|. The Lienard-Wiechart po- 
tentials are employed in Sect. (jIII|) to take into account 
the retardation effects. Sect. pil|) also includes a brief 
discussion to illustrate conditions under which the differ- 
ences between the old (wrong) solution and the new (cor- 
rect) solution are important. Discussion of the correspon- 
dence to analagous gravitational problems is included in 
Section pV[) . A conclusion is included to summarize the 
main findings and to introduce other problems in elec- 
tromagnetism where such effects could be important. 



II. THE CLASSIC IMAGE PROBLEM 

The classic image problem in electromagnetism in- 
volves a grounded infinite conducting plane and a point 
charge. The force of attraction between the point charge 
and the plane is identified using an equivalent picture of 
an image charge equidistant behind the plane of opposite 
sign. The field configuration in the half-space of the point 
charge is mathematically equivalent to the picture where 
an image charge of equal magnitude and opposite sign 
resides equidistant behind the plane. Hence, the field 
configuration of the point charge and the plane is math- 
ematically equivalent to the the field between the point 
charge and the image charge in the space of interest. 

If released from rest, the charge begins to move toward 
the plane, which also moves the image charge closer to 
the midpoint between the charge and the image. Conse- 
quently, the force of attraction between the charge and 
the image rapidly increases. The expectation in the text- 
books is that Newton's 2nd law can be employed to write 
a second-order differential equation to be solved for f(t). 
Typically, one is asked to show that the time to reach the 
plane isj!] 
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wher e q is the charge, m is the mass and d is the initial 
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permittivity of free space appears as ep which implies a 
particular system of units (SI) . 

Start with a point charge q of mass m released from 
rest a distance d from an infinite grounded conducting 
plane (sec FigH]). An image charge q' = -q placed a dis- 
tance d behind the plane will result in the correct bound- 
ary condition for the electric potential V{0) = and 
the correct field configuration for a; > 0. Consequently, 
by uniqueness theorem arguments, one knows that the 
image configuration is mathematically equivalent to the 
original problem. As soon as the original charge moves, 
so too will the image charge (very much like two gravitat- 
ing masses approaching one another; however, the image 
charge is not a physical charge so it does not introduce 
kinetic energy into the picture). 



q, m 



FIG. 1. Problem setup: Point charge released from rest a 
distance d from a grounded infinite conducting plane. 



The one-dimensional equation of motion under these 
conditions is (positive x to the right) 



dr- 



(2) 



where the initial conditions in time are a;(0) = d and 
w(0) = 0, with the endpoint condition x{t) ~ 0. Upon 
integrating once and applying the initial condition(s) one 
can write 
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dx 
'dt 



(3) 



Already here there is a problem! When the charge 
reaches the plane the value of a: becomes 0. Therefore, the 
particle's speed becomes unbounded — clearly unphysical. 

Still, if one proceeds to integrate once again to find an 
expression relating time t and position x, one finds 



t = T 



1 sin 



(4) 



Setting a; = identifies the time to reach the plane, but 
one must remember that particle's speed has become in- 
finite at that moment. 

How can this be corrected? The first step toward a 
complete solution is to use the relativistic expression for 
the particle's momentum p = mjv. Then, changing the 



time derivative to a spatial derivative through dt — vdx, 
we can write 

-g2 d{jv) 



Wneox^ 
—q^ dx 



dt ' 
v"f'^ dv . 



(5) 



The Lorentz factor 7 = (1— w^/c^)"^/^ appears charac- 
teristically to modify the dynamics appropriately when 
the kinetic energy is not small compared with the rest 
energy. This expression can be integrated to find 



(l-g)V2fa(l-f) 
l{l-^)+x 



(6) 



where a characteristic length £ = g^/(167reomc^) ap- 
pears. Notice that v{d) = and v{{)) = — c (with the 
minus sign pointing to the direction) as it must to avoid 
superluminal transport. For a specific choice of charge, 
mass, and initial separation [c?/(cr) = 0.5, where d/{cT) 
= {l/Tr)^y8£/d], a comparison is presented in Fig. [2] of 
the relativistic versus nonrelativistic expressions for the 
speed. For this particular configuration, the difference 
seems significant. 
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1 -x/d 

FIG. 2. Magnitude of velocity divided by c versus distance 
traveled (that is, 1 minus the position divided by d) for a 
point charge released from rest near a grounded conducting 
plane. The solid curve (Eq. [6]) describes relativistic dynamics 
and the dashed curve corresponds to nonrelativistic dynamics 
(Eq.[3l). 



The next task is clearly to integrate Eq. ([6]) to find an 
expression for i as a function of x. Specifically, we must 
integrate 



ct 



■.dx'. 



(7) 



/^2(i_^) +2fe'(l-^) 
The integral can be done in closed form and written as 
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A position-time analysis is now possible by comparing 
Eqs. (H]) and ([5]). The comparison is shown in Fig. 
with the solid curve representing relativistic dynamics 
and the dashed curve displaying the nonrelativistic ex- 
pression. The same choice is made here for the di- 
mcnsionless parameter as was made previously [namely, 
d/{cT) ^ 0.5]. 




t/x [d/(CT)] 

FIG. 3. Position (distance traveled in units of d is "1— x/rf") 
versus time for a point charge released from rest near a 
grounded conducting plane. The dimensionless parameter is 
chosen to be d/(cr) = 0.5. The solid curve corresponds to rel- 
ativistic dynamics and dashed curve uses the nonrelativistic 
expression. 



The slope of the tangent line is of course a measure of 
the magnitude of the instantaneous velocity. However, 
it cannot be read off directly since there is a factor of 
d/{Tc) in the ratio. And yet, notice that the dashed 
curve exhibits a slope for its tangent line near x = 
that is very large. Numerically the magnitude of the 
velocity exceeds 1 (in units of c) for t/r > 0.62 in the 
dashed curve. At t/r = 0.8 the speed exceeds 2, and 
it exceeds 4 at f/r of 0.96. Furthermore, we point out 
that the rclativistically consistent curve displays a slope 
corresponding to a speed less than 1 throughout, only 
does the speed approach 1 at the point x = Q. 



Now that we have the expression for i as a function of 
a;, it is a straightforward exercise to take x to zero and 
identify the time to reach the plane. By setting a; — > 0, 
then t ^ and we have, 




(9) 



where r is the original nonrelativistic expression given by 
Eq. ([T|). We remark that if c — > oo, the right hand side 
approaches 1, so T— >■ r. 

A comparison between the nonrelativistic expression 
for the time to reach the plane as compared with the fully 
relativistic expression is ultimately useful to be able to 
identify the conditions under which the original textbook 
expression fails. This comparison is deferred until after 
a discussion of the field-retardation effects. 



III. FIELD RETARDATION EFFECTS 

The discussion of the previous section improves upon 
the nonrelativistic analysis which is typically used to find 
the time to reach the plane. However, it too is lacking in 
some respects. When the charge is a distance x from the 
plane and moving toward the plane, the force it experi- 
ences is different from a static Coulomb picture because 
the "electromagnetic news" travels at the speed of light. 
Hence, the electric field it experiences must be related to 
the electric and magnetic potentials at some earlier time 
(the so-called retarded time). Simply put, it takes a time 
x/c for information on the particle dynamics to reach the 
plane. The charge configuration on the plane is updated 
(which means the image charge location is updated) and 
then that information requires a time 2:/c to get back 
to the point charge. This type of physical circumstance 
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is describable by the Lienard-Wiechart potentials for a 
moving point charge [1, 

Suppose the charge is located at position x and moving 
with a velocity —v (again, positive x to the right). The 
general problem to calculate the retarded potentials is 
nontrivial, but when the motion is one-dimensional it 
simplifies significantly. In particular, the electric field 
experienced by the charge in this case is reduced by two 
powers of 7. It can therefore be written as 



10' 
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E = 



(10) 



The force drawing the charge toward the plane is there- 
fore velocity dependent and it is reduced as compared 
with the static case. Consequently, the rclativistic and 
properly time-delayed equation of motion is 



-q 



V 



d{'-fv) 
' dt 



167reox2 J dt ' ^^^^ 

Remarkably, this too can be integrated to find the veloc 
ity as a function of position. It becomes 



dx 



xd 



xd + U{d - . 



2/3 



(12) 



Another integral is required to identify a relationship be- 
tween position and time in this relativistically consistent 
and properly casual description of the dynamics. Unfor- 
tunately, a closed form expression has not been found. 
Instead, the integral will be carried out numerically for 
the case of interest here. That is, the time T to reach 
the plane will be computed and displayed graphically to 
help clarify the size of the effects. In particular, the time 
T is 



T 

T 



dz 



.(13) 



1 



Time T is normalized here by the original nonrelativistic 
time. A ratio of 1 indicates that the original nonrelativis- 
tic expression is completely satisfactory. Deviations from 
1 indicate that relativistic effects and/or field retardation 
effects become important. 

The time to reach the plane is again a function of the 
dimensionless parameter d/(rc). In Fig. (|4]) the full so- 
lution is plotted as a function of d/{Tc). One must keep 
in mind that the dimensionless parameter is related to 
the charge, mass, and initial separation through 

- = -\/? - -^^=- (14) 
Tc TT y d TTy'2TTeomc^d 

Graphical results indicate that so long as the dimen- 
sionless parameter d/{Tc) is less than roughly 0.1, then 
the nonrelativistic expression for the time is completely 
adequate. For values of the parameter greater than 1 
the relativistic and retardation effects become very im- 
portant. What remains to be considered is to identify 
the physical situations that correspond to each of these 
regimes. 
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FIG. 4. Time to reach the plane versus scaled initial sepa- 
ration. Solid curve includes relativistic corrections, dashed 
curve includes relativistic correction plus retarded (Lienard- 
Wiechart) potentials. 



IV. DISCUSSION 

There is of course some concern about using a solu- 
tion for T that involves super luminal transport. How- 
ever, we have identified the crossover point where the 
relativistic expression must be used (see Fig IH) . If a par- 
ticular application falls to the left of the crossover point 
on Fig. (j4]), then the nonrelativistic formulation is com- 
pletely adequate. It will be useful to consider a few ex- 
ample configurations to get a feel for this dimensionless 
parameter d/{Tc). Typical laboratory-type charges are 
microcoulombs, with typical masses in the milligram-to- 
gram range, and initial separations in the neighborhood 
of tens of centimeters. 

For the typical configuration described above, the im- 
portant parameter comes out to be of the order 10~^ 
as can be seen from the information in Table HI Conse- 
quently, the nonrelativistic solution is completely suffi- 
cient for this case. We emphasize that even though the 
nonrelativistic expression for the total time to reach the 
plane is valid, it does not mean the particle was com- 
pletely nonrelativistic. Rather, the time spent moving 
relativistically was such a small part of the total time 
that relativistic corrections are unimportant. It is espe- 
cially reassuring to notice that the number 10^^ is several 
orders of magnitude smaller than the region in Fig. ^ 
where the relativistic and retardation effects begin to be 
important. 

Next we explore the case of an electron released from 
rest a distance 10 cm from a grounded conducting plane. 
Again from the table we extract the important dimen- 
sionless parameter for this case to be a number of the or- 
der of 10~® (even smaller than the previous case!). This 
too is completely legitimately described by nonrelativis- 
tic formulae. The argument that saves the nonrelativis- 
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q 


d 


m 


d/{rc) 


ImC 

electron 
electron 
proton 


1 (cm) 
0.1 (meter) 
1 (fm) 
10 (fm) 


1 milligram 

me 
me 
nip 


~ 10"'^ 
7.5x10"^ 
~ 1 

~ 6x10"^ 



TABLE I. Comparison of different physical configurations 
and the resulting dimensionless parameter d/(rc) which de- 
termines the region of applicability of relativity and field- 
retardation effects. 

tic expression is that even though the last tiny fraction 
of the trip to the grounded plane is superluminal, the 
charge spends so little time moving faster than light that 
the the deviation from the "correct" time is immeasur- 
ably small. So, once again, the nonrelativistic expression 
is sufficient. 

Lastly, we investigate extreme cases where relativistic 
and retardation effects could be important. Consider an 
electron released from rest a distance 1 femtometer from 
a grounded conducting plane. The important parameter 
determining the dynamics is again the number d/{Tc) 
which comes out to be of the order 1 for this configura- 
tion. The nonrelativistic expression will fail noticeably 
for this case. Not only is relativity important for the 
dynamics, but field retardation effects are very impor- 
tant. A proton released 10 fm from the plane results in 
a parameter d/{Tc) of 10^'^. 

Before concluding, it is useful to extend the solution 
presented above for the electrical charges to the case of 
two gravitating point masses. Consider two point masses 
m initially separated by a distance 2d. We have immedi- 
ately solved that problem as well using a correspondence 
between charge and mass, and between the inverse of the 
permittivity of free space and Newton's universal gravita- 
tional constant. The above dynamical expressions relat- 
ing position, velocity and time hold for the case of gravity 
while the relevant dimensionless parameter becomes 

d ^/2 , , 

TC TT y dC' 

where G is Newton's universal gravitational constant. 
For most configurations this parameter is extremely 
small. Hence, once again, the nonrelativistic expression 
for the time to reach the center of mass is completely safe 
to use. 

V. CONCLUSION 

The classic image problem of a point charge near a 
grounded conducting plane invites the notion of a math- 
ematical tool of replacing the grounded plane with a fic- 
titious charge (a so-called image charge) which, together 
with the original charge, produces the same electric po- 
tential in the region of interest and respects the boundary 
conditions. The existence and uniqueness theorems indi- 
cate that the scalar potential and therefore the (vector) 



electric fields must be equivalent. It is a beautiful solu- 
tion to an otherwise difficult problem. 

The follow-up problem involves computing the time it 
takes the charge to reach the plane when the charge is re- 
leased from rest. A straightforward (although incorrect) 
analysis from nonrelativistic classical mechanics results 
in an expression for the time which is proportional to 
the initial separation to the three-halves power, also ro- 
portional to the square root of the mass, and inversely 
proportional to the charge. We have pointed out in this 
article the serious formal flaw in that proposed solution. 
The flaw is linked to the observation that the velocity of 
the charged particle exceeds the speed of light for a small 
portion of its trip toward the grounded plane. It is very 
important to sort out when the solution is still valid and 
acceptable, and when it is noticeably inadequate. This 
has been the purpose of our article. 

We have shown that the relativistic version of the prob- 
lem is completely tractable and results in subluminal 
dynamics which changes the expression for the time to 
reach the plane by monotonically increasing the result as 
a function of the dimensionless parameter d/irc). 

A more sophisticated approach to point charge dynam- 
ics near a grounded conducting plane must include field- 
retardation effects. This has also been added to the so- 
lution discussed in this paper. Using time-delayed fields 
(Lienard-Wiechart potentials) increases further the time 
required for the particle to reach the plane. Both effects, 
relativistic and retardation influences, are easily under- 
stood in terms of their tendency to increase the time for 
the particle to reach the plane. 

Finally, we briefly mention other problems where these 
same effects could be important. First, a point charge 
q of mass m near a fixed line charge —A experiences 
a force that is proportional to l/r. The nonrelativistic 
time to reach the line charge can be computed asQ r = 
■ndyJeQm/{\q). However, the charge goes superluminal 
as it reaches the line charge. A point dipole p of mass 
m near a fixed line charge —A experiences a 1/r^ force. 
Using nonrelativistic dynamics again we find the time to 
reach the line charge is t = ij: d / A) ireQmd / {pX) . Fur- 
thermore, a point charge q of mass m near a point dipole 
p of mass m will experience an inverse cube force (1/r^). 
The nonrelativistic time for the two to meet at the mid- 
point (each of which is initially a distance d away from 
the midpoint) is r = d'^ yjlixt^mj (j)q ). And lastly, two 
identical point dipoles p with mass m will experience 
a force proportional to l/r*. The nonrelativistic expres- 
sion for the time to reach the midpoint isQ r = (2/3)K(i^ 
y/neonid /p, where the constant k can be computed as 



.7r/2 
K= / 








Jo 





In each case above, analytic or numerical expressions for 
velocity as a function of position and position as a func- 
tion of time can be identified. Relativistic corrections 
and field retardation effects are equally readily obtain- 
able. For increasingly steeper r-dependence, one might 
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expect that the effects could be increasingly more impor- 
tant. The small-distance behavior determines the relative 
importance of the effects, and hence the dipole-dipole is 
the most sensitive to these effects. 



ACKNOWLEDGMENTS 

This work was supported in part by the National Sci- 
ence Foundation under grant No. PHY-0555521. 



[1] R. Resnick, D. Halliday and K.S. Krane, Physics, Vol. 1, 
Fifth Edition, J. Wiley & Sons, Inc., 328 (2002). 

[2] S.T. Thornton and J.B. Marion, Classical Dynamics of 
Particles and Systems, Fifth Ediction, Brooks/Cole, 205 
(2004). 

[3] B.F. Bayman and M. Hammermesh, A Review of Under- 
graduate Physics, J. Wiley & Sons, Inc., 117 (1986). 

[4] D.J. Griffiths, Introduction to Electrodynamics, Third Edi- 
tion, Prentice Hall Inc., 155 (1999). 

[5] S. Palit, Principles of Electricity and Magnetism, Alpha 
Science International Ltd., 48 (2005). 

[6] J.D. Jackson, Classical Electrodynamics, Third Edition, 
661 (1999). 

[7] D.J. Griffiths, Introduction to Electrodynamics, Third Edi- 
tion, Prentice Hall Inc., 435 (1999). 

[8] The expression to find r is proportional to / = 
/o dV\/V ■ Square the integral, which becomes a dou- 



ble integral over rj and ly. Change variables to u = rj + u 
and V — rj — v. The area element drjdu becomes 2dudv. 
The resulting integral is soluble with the result coming out 
as / = ^/7r/2- 

[9] In this case, r is proportional to k = 
dr] e"^^''^'''/Vl - e^') . Expanding and integrat- 
ing term by term gives 



"Z^ (6n- l)(n- l)!2"-i 
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The same result comes from numerically integrating 
Eq. (HSl). 



